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Abstract 

We derive a general set of Poor Man's scaling equations and analyze the stability of 
the Luttinger state in a system composed of a finite number N of one dimensional spinless 
fermionic chains, coupled through a general two body interaction. The effect of processes 
with momentum transfer parallel to the Fermi surface in destroying massless states is inves- 
tigated. It will be shown that there are two processes competing: one in which two electrons 
exchange chains and the other in which they jump into a same chain. When periodic bound- 
ary conditions in the transverse direction are taken into account this competition leads always 
to massive states (except in hyperplanes of the phase diagram) , a well known example being 
the generalized sine-Gordon model. If instead open boundary conditions are taken, mass- 
less states are possible but due to this competition the system is placed near instabilities. 
We argue that this kind of analysis has relevance for understanding the instabilities of 2D 
fermionic systems. 
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The many unexplained features of high-Tc superconductors have motivated an intense study 
of low dimensional interacting fermions(l|. A renormalization group (RG) approach then de- 
veloped allowed a reformulation of the Fermi liquid concept [Q]. However this analysis seems to 
rely heavily on the isotropic nature of the Fermi surface(FS). Indeed, for a two-dimensional(2D) 
circular FS one can neglect, to some extent, the effect of scattering processes with momentum 
transfer tangential to the Fermi surface. This is because the FS topology highly restricts con- 
tributions to the zero sound channel by a factor of (A/hp) 2 (where A is the momentum cutoff 
used in the RG procedure), whereas the Cooper channel receives contributions proportional to 
A/fcp. Thus we can neglect the zero sound contributions to the one loop renormalization of the 
four fermion vertex and then we are left with a simple differential equation that diverges only if 
attractive interactions are present, indicating the appearance of a superconducting instability. 
However for a FS with large flat portions (nested), like in the case of a square FS, this picture 
may not apply. Then the scattering processes with momentum transfers parallel to the FS and 
involving electrons on parallel sides of the FS, give supplementary contributions and may play 
an important role in destabilizing massless regimes. To understand better the effect of these 
contributions, we suggest, in a first reasonable approximation, the study of a model of ID chains 
coupled through a general two body interaction. 
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It is worth noting that this kind of model has quite a large range of applications, for it can 
be related (even if approximately) to ID spin chains through the Jordan- Wigner transformation 
and to systems of coupled Hubbard chains^]. Even though, these are particular cases of our 
model since for these models the chains are coupled through density-density interactions, and 
the densities are operators defined on each chain. Consequently the number of electrons on 
each chain is conserved. By contrast, this is not a requirement for our model, which has this 
new degree of complexity. In particular, some commonly used techniques, like ID bosonization, 
find inconsistencies if one takes account of the fermionic nature of the field operators correctly, 
through the Haldane ladder operators g. For this reason we derive a general set of RG 
differential equations based on the Poor Man's Scaling (PMS) method, thus avoiding the use of 
the bosonization tools whenever they lead to incongruities. 

Consider an array of N spinless fermionic chains, aligned along x. As in the ID case we 
separate the field operator in slow and fast components, ^(x) = l/\/~NJ2p=± 2~2 n =l ex P(^pk n • 
x )V , p,n( x ) ) where pk n is the momentum Fermi surface point for the chain n and branch p. We 
take all (k n ) x equal to kp, a linearized dispersion relation near the FS, and neglect any Fermi 
velocity modulation. A general density-density interaction can be written as: 

1 S 6 

W C = ^EEE G s( i ~ 2' J + 2^ R i+8/2,k+q L t-5/2,k'-g L j+6/2,k> R%-S/2,k i 1 ) 
i,j,8 q k,k' 

where here L~^ k is the creation operator for a left mover on chain j and with momentum k along 
the chain. If we take periodic boundary conditions (PBC) along y, then the couplings depend 
only on j — i = A; Gs(I, J) — * Gg(A), with J = j + 6/2,1 = i — 5/2. From (JTJ) we can see 
that due to the inversion symmetry (i<->j), hermiticity and boundary conditions, these couplings 
must verify the conditions: 

G s (I,J) = G_ s (J,I)=G_ s (I + S,J-5) = G„ s (N + l-I,N+l- J) (no PBC) 
Gs(A) = G^(A) = G 5 (-A) = G^(-A) = Gs±n(A) = G 5 (A±N) (PBC) [Z> 

Note that the chain index can also be regarded as a spin index. In this way we can establish 
a dictionary between the usual g-ology couplings |?]] and the G's, for a model with 2 coupled 
chains: 

G (1)/2 = S2± Gi(0)/2 = -su. G (0)/2 = 0a|| -Sin Gi(l)/2 = (3) 

Notice that according to the formulation (|l]), the backscattering processes are now seen as 
forward scatterings with spin-flips, and thus within each chain there are no backscatterings (see 
figure la ). This formulation will allow us to obtain a compact set of RG equations, which is 
particularly useful in treating the N-chain problem. 

The PMS equations are derived by requiring invariance of the vertex calculated in a one 
loop expansion. There are only two diagrams giving logarithmic contributions, both with the 
same magnitude (po/2 [ln(u)/D) — in/2] ) but having opposite signs. Here po is the density of 
states for one spin direction. For convenience we will redefine from now on the couplings as 
the adimensional quantities Gg(A) — > Gg(A)/irVF- The RG flow equations for models with and 
without PBC are then: 

= [G a (A + 5-a)G 5 - a (A- a )-Gs- a (A)G a (A)} (PBC) (4) 

a=0 
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Q/~f (J J\ -1 N — 1 

= 2N ^ [G a (I,J)G S - a (I + a,J-a)-G S -a(I + a,J)G a (I,J- 5 + a)](5) 

a=—N+l 

where in (^) the boundary conditions impose that only couplings with 1,1 + 5, J and J — 5 
between 1 and N are non-zero. The reader may verify that all the RG equations found in ^ 
are reproduced by @, via the dictionary (||). The generalized flow invariant is ^a^o^)' It 
is important to call attention to the compactness of the formula (|3]). Note that if one takes 5 
and A as vectors instead of scalars, then the couplings can be properly chosen to obtain the RG 
equations for models of coupled chains of electrons with spin. 

Let us analyze first systems with PBC. For N=2, bosonisation of the model 
gives[||: 7i = Htl+ ~Hgsg, where Htl is the Tomonaga-Luttinger hamiltonian, and Hgsg = 
l/(2vra) 2 / dx -Gi(0) cos (\/80-) +Gi(l) cos (V86_ ) . The model consists of the sine- 
Gordon model, plus a cosine term in the dual field. For this reason it is called the generalized 
sine-Gordon model (GSG). The GSG model is always massive because it involves dual operators, 
which thus have inverse scaling dimensions. There can only be gapless states in hyperplanes 
of the phase diagram. A generalization to the N-chain problem leads to the same conclusion. 
A simple way of seeing this uses a linearization of the equation (0), around the possible fixed 
points. Writing Gg(A) = G|(A) + 55(A) and linearizing around the Luttinger Liquid fixed 
point, we obtain simply: 

dg s (A) = 1/2JV (G* (A + S) + G* (A -5)- 2G* (A)) g s (A) (6) 

We realize that we can always have two symmetrical eigenvalues if N is even: for A = 5 = N/2, 
A = 1/2N [G (0) - G (N/2)} whereas for A = 0, 5 = N/2, A = 1/2N [G (N/2) - G (0)]. This 
corresponds to a competition between two quite different kinds of processes: in the first case two 
electrons jump within the same chain, whereas in the second they exchange chains. Thus we 
conclude that this competition has a dramatic effect in destabilizing gapless states in a system 
where PBC are taken. An example of systems behaving in this way are the systems of chains 
coupled by transverse hopping terms p. 

Next we investigate the case where no PBC are considered which is the model with more 
practical relevance. Here the boundary conditions do not allow processes like G_/v/ 2 (N/2) and so 
we may think that this kind of competition is not present. We show next that this competition 
is indeed softened, but that it is still present although in a slightly different fashion. Consider 
a system of four coupled chains (no PBC). A simple argument states that if the only non- 
zero couplings are Gi (2, 2), Gi (2, 3), Go (2, 2) and Go (2, 3) (and all the other couplings that by 
symmetry are equal to these), then apart from multiplicative factors, the RG equations have the 
same structure as those of the GSG model. Intuitively this arises because we can look at the first 
chain as the image of the third and the fourth as the image of the second, like if we had a model 
with PBC. Thus we showed that in this region of the phase diagram a four chain model with 
open boundary conditions has the same properties as the GSG model. This subspace is stable 
and doesn't really tell us much about the effect of the omitted couplings. However this provides 
us with a picture for a kind of competition present between couplings of the form Gs (1, 1 + 5) 
and G s {I, I) • 

Consider a model with three spinless chains. We rename the couplings as shown in figure 
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lb). Then the RG equations look like: (the factor 1/N is omitted ) 

d(A + C + D) = d(B + 2C) = d{B - 2C) = 4(£ 2 - F 2 ) 

d{D + A) = (E 2 - F 2 ) dE = E(B -2C + G + D) (7) 

8G = E 2 - F 2 + 2G{D -A) dF = —F(B -2C + G + A) 

Linearizing around the various fixed points we can verify that the only stable fixed points are 
of Luttinger type (E=F=G=0) and lie in the region A < 2C — B < D. This seems to hold 
true for an arbitrary number of chains as the number of conditions for the eigenvalues of @ 
is proportional to N 3 (the order of the number of different Gs^o(I, «/)), whereas the number of 
relations we can have among the G$=o(I,J) couplings goes as N 4 . Thus finding it possible to 
have a weak coupling theory for N=2 and 3, it looks likely that a weak coupling theory exists 
for any N. 

Another point worth noticing is how the E and F couplings do infact compete. They flow 
independently of their sign and almost in opposite directions. On the other hand, the 5 = 
processes depend only on the relative magnitude between the E and F processes, which is a clear 
manifestation of their competition. For positive bare couplings the system can reach two kinds 
of strong coupling regimes. If F^ +oo and E — > 0, then A,B,G^ +oo, D,C^ — oo (Regime I). 
If E^ +oo and F -> 0, then G,D,B^ -oo, A,C^ +oo (Regime II). 

We are interested in knowing whether what we found for a system of 3 chains provides us 
with a picture for systems with a higher number of chains. Here we argue that this can indeed be 
the case for a large class of systems. We look at systems close to a generalized Luttinger liquid 
state and for strictly repulsive interactions. We want to study the influence of small momenta 
transfer processes (tangential to the FS) in destabilizing the massless Luttinger state. For this 
reason we analyze bare coupling functions of the general form Gs(I, J) = g(\5\) f(J — 1 — 5). 
The function g controls the dependence of the couplings on the momentum transfer. We take 
it as a positive monotonic (decreasing) function. The forward scattering models used in studies 
of 2D bosonization correspond to the case where g(\5\) is a delta function at 5 = 0. It's easy to 
show that G respects all the necessary symmetries, provided / is an even function. Here we will 
present some numerical results. An analytical study of the equation (l|) for this class of bare 
coupling functions can also be achieved but this is left for future publication. 

In the following numerical results we considered systems of five chains. We take <7(|<5|) = 
exp (— 5 2 j2a) . The form of g doesn't seem to play an important role so long as it is even 
and peaked only at 5 = 0. The control parameter a may be seen as anisotropy dependent: 
we can expect a curved FS to favor small values of a. For the function / we take / (x) = 
ao (l + &o ( X /N) 2 ^ , like in a Taylor expansion. The coefficients ao and &o are chosen so that / 
remains positive, as we are interested in studying the effect of repulsive interactions. In principle 
it may look physically more reasonable to take / as a decreasing function. However one should 
also remember that / is already a renormalized function where other effects like those coming 
from phonons and anisotropy are included. 

The results can be summarized as follows. When / is monotonous and decreasing then 
no matter how small a is, we fall in a regime where all chain exchange processes diverge to 
+oo and the jump within the same chain processes decrease (figure 2a ). This corresponds to a 
generalization of the regime I (F— » +oo), in the three chains model. The competition between 
the two classes of couplings is obvious. Then we also clearly observe that all Go (n, m) , with 
n ^ m, diverge to — oo whereas all the diagonal elements Go (n,n) , diverge to +oo. When / 
is monotonous and increasing we can have several regimes. Like in the three chain model a 
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Luttinger state is possible for small enough a. By increasing a a massive regime appears similar 
to regime II (E — * +00) of the three chain model (figure 2b ). Now some processes with jumps 
within the same chain may also diverge to —00. The remaining couplings behave in the same 
way as in the three chain system. The competition between each pair of couplings G5 (I, 1 + 5) 
and G$ (I, I) , is still present but works in different ways and magnitudes for each / and 5. If a 
is further increased then we recover regime I, where the jumps within the same chain processes 
strongly decrease (figure 2c ). In the three chain model it is easy to see why this happens: by 
increasing a we increase the coupling G, which is a chain exchange process. Thus this class of 
processes is favoured and manages to overcome the jump within the same chain processes. This 
analysis seems to apply for any finite N. We would like to stress that the three chain model 
seems to reproduce much of the physics of a more general N chain model. 

Finally, we consider a 2D electron gas, where flat portions are introduced in a circular 
FS (figure 3). We take N patches in the nested regions, and M in the curved portions, with 
N/M^; 1. Due to the geometry of the FS, it is more appropriate to parametrize the var- 
ious patches in terms of their angular positions. This is achieved through the substitution 
G$ (I, J) — ► G$ (I, N + M + 1 — J). The symmetries on the square restrict the number of in- 
dependent couplings: 0< 5 < N + M, 1< I,J< N + M. Using Shankar's arguments Q, we 
neglect any zero sound contributions to a vertex involving electrons on the curved portions. 
This corresponds to introducing a factor 81^61+5^ on the second term of (|5|) (<fi stands for 
the nested region). Also, the only couplings with electrons on these regions have the BCS form 
Gs(I,I)- With these modifications, the equation (||) is still valid, if we regard the couplings 
as already renormalized by the finite 2D density of states. We linearized the RG equations 
around the Luttinger fixed point for a system with N=3, M=0 and concluded that the Luttinger 
regime remains stable if a further condition on the couplings B,D>0, is verified. The other 
strong coupling regimes do also exist. If the couplings involving scatterings not restricted to 
a nested region are sufficiently strong, then the regime II leads to a new regime, where the E 
coupling goes to -00 (figure 3). This state seems to be a good candidate for a superconductive 
state. In fact, in this regime F— > 0, so that the main corrections to the vertices come from the 
Cooper channel, and are due to effective attractive interactions. This will be studied in detail 
elsewhere 

We should point that the possibility of generating ordered phases due to nesting is not 
a new concept: it remains the main explanation for itinerant antiferromagnetism,it has been 
suggested to explain the HTCSC, and it may explain the stability of some CDW phases recently 
observed ||]. 

We acknowledge discussions with A.M.Tsvelik. FVA thanks discussions with R. Chitra, L. 
Laloux, B. Delamotte, J.M.Robin and financial support from EEC under the grant ERBICH- 
BICT941242. 
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Figures Captions 




Figure 1. a) One-dimensional interacting electrons with spin can be seen as two spinless chains 
in interaction. The backscattering process (dashed arrows), where two electrons change from one 
side of the FS to the other, is now seen as a spin- flip process (solid arrows). In b) all scattering 
processes involved in the three chain model with no PBC are shown. 
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Figure 2. Strong coupling regimes: a) f (x) = 0.4 ^1 — 1/2 (x/N) 2 ^j and a = 0.5; 
b) / (x) = 0.4 (l + (x/iV) 2 ) and a = 0.5; c) same / (x) but <r = 2.5 
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Figure 3. For a system with N=3, M=14, bare couplings A=B=C=D=0.1, E=F=0.08, 
G=0.02 and the remaining couplings equal to 0.12, the Fourier modes V/ of J2i (I, I) diverge 
to -oo, whereas F^ 0. 
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